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Abstract 

A canonical result about satisfiability theory is that the 2-SAT problem can be solved in 
linear time, despite the NP-hardness of the 3-SAT problem. In the quantum 2-SAT problem, 
we are given a family of 2-qubit projectors n^- on a system of n qubits, and the task is to 
decide whether the Hamiltonian H = ^ 0-eigenvalue, or it is larger than l/n“ for 

some a = 0(1). The problem is not only a natural extension of the classical 2-SAT problem to 
the quantum case, but is also equivalent to the problem of finding the ground state of 2-local 
frustration-free Hamiltonians of spin a well-studied model believed to capture certain key 
properties in modern condensed matter physics. While Bravyi has shown that the quantum 
2-SAT problem has a classical polynomial-time algorithm, the running time of his algorithm is 
O(n^). In this paper we give a classical algorithm with linear running time in the number of 
local projectors, therefore achieving the best possible complexity. 
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1 Introduction 


Various formulations of the satisfiability problem of Boolean formulae arguably constitute the center 
piece of classical complexity theory. In particular, a great amount of attention has been paid to 
the SAT problem, in which we are given a formula in the form of a conjunction of clauses, where 
each clause is a disjunction of literals (variables or negated variables), and the task is to find a 
satisfying assignment if there is one, or prove that none exists when the formula is unsatishable. In 
the case of the A;-SAT problem, where A; is a positive integer, in each clause the number of literals 
is at most k. While /s-SAT is an NP-complete problem [Coo71, Kar72, Lev73] when k > 3, the 
problem 2-SAT is well-known to be efficiently solvable. 

Polynomial time algorithms for 2-SAT come in various flavors. Let us suppose that the input 
formula has n variables and m clauses. The algorithm of Krom [Kro67] based on the resolution 
principle and on transitive closure computation decides if the formula is satishable in time 0{rfi) and 
finds a satisfying assignment in time 0{n^). The limited backtracking technique of Even, Itai and 
Shamir [EIS76] has linear time complexity in m, as well as the elegant procedure of Aspvall, Plass 
and Tarjan [APT79] based on computing strongly connected components in a graph. A particularly 
simple randomized procedure of complexity O(n^) is described by Papadimitriou [Pap91]. 

For our purposes the Davis-Putnam procedure [DP60] is of singular importance. This is a 
resolution-principle based general SAT solving algorithm, which with its refinement due to Davis, 
Putnam, Logemann and Loveland [DLL62], forms even today the basis for the most efficient SAT 
solvers. While on general SAT instances it works in exponential time, on 2-SAT formulae it is of 
polynomial complexity. 

The high level description of the procedure for 2-SAT is relatively simple. Let us suppose that 
our formula (p contains only clauses with two literals. Pick an arbitrary unassigned variable Xi and 
assign Xi = 0. The formula is simplified: a clause {xi V Xj) becomes true and therefore can be 
removed, and a clause {xi V Xj) forces Xj = 1. This can be, in turn, propagated to other clauses to 
further simplify the formula until a contradiction is found or no more propagation is possible. If no 
contradiction is found and the propagation stops with the simplified formula (po, then we recurse 
on the satisfiabilty of (pQ. Otherwise, when a contradiction is found, that is at some point the 
propagation assigns two different values to the same variable, we reverse the choice made for Xi, 
and propagate the new choice Xi = 1. If this also leads to contradiction we declare cp unsatishable, 
otherwise we recurse on the result of this propagation, the simplihed formula (pi. 

There is a deep and profound link between fc-SAT formulas and /c-local Hamiltonians, the central 
objects of condensed matter physics. A A;-local Hamiltonian on n qubits is a Hermitian operator 
of the form H = Yll^i where each hi is by itself a Hermitian operator acting non-trivially on 
at most k qubits. Local Hamiltonians model the local interactions between quantum spins. Of 
central importance is the minimal eigenstate of the Hamiltonian, known as the ground state, and 
its associated eigenvalue, known as the ground energy. The ground state governs much of the 
low temperature physics of the system, such as quantum phase transitions and collective quantum 
phenomena [Sac07, VLRK03]. Finding the ground state of a local Hamiltonian shares important 
similarities with the /c-SAT problem: in both problems we are trying to hnd a global minimum 
of a set of local constraints. This connection with complexity theory is of physical significance. 
Indeed, with the advent of quantum information theory and quantum complexity theories, it has 
become clear that the complexity of finding the ground state and its energy is intimately related to 
its entanglement structure. In recent years, much attention has been devoted into understanding 
this structure, revealing a rich an intricate behaviour such as area laws [ECPIO] and topological 
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order [Kit03]. 

The connection between classical fe-SAT and quantum local Hamiltonian was formalized by 
Kitaev [KSV02] who introduced the A:-local Hamiltonian problem: one is given a fe-local Hamiltonian 
H, along with two constants a < b such that b — a > 1/n" for some constant a. It is promised that 
the ground energy of H is at most a (the YES case) or is at least b (the NO case), and the task 
is to decide which case holds. Broadly speaking, given a quantum state {ijj), the energy of a local 
term {il;\hi\ip) is a measure of how much \ip) “violates” hi, hence the ground energy is the quantum 
analog of the minimal number of violations in a classical /c-SAT. Therefore, in spirit, the fc-local 
Hamiltonian problem corresponds to MAX-A:-SAT, and indeed Kitaev has shown [KSV02] that 5- 
local Hamiltonian is QMA-complete, where the complexity class QMA is the quantum analogue of 
classical class MA, the probabilistic version of NP. 

The problem quantum A:-SAT, the quantum analogue of fc-SAT, is a close relative of the fc-local 
Hamiltonian problem. Here we are given a A:-local Hamiltonian that is made of /c-local projectors, 
H = Qij asked whether the ground energy is 0 or it is larger than b = 1/n" for 

some constant a. Notice that in the YES case, the energy of all projectors at the ground state is 
necessarily 0, since by definition, projectors are non-negative operators. Classically, this corresponds 
to a perfectly satisfiable formula. Physically, this is an example of a frustration-free Hamiltonian, in 
which the global ground state is also a ground state of every local term. Bravyi [Brail] has shown 
that quantum fc-SAT was QMAj^-complete for k > 4, where QMA^ stands for QMA with one-sided 
error (that is on YES instances the verifier accepts with probability 1). The QMA;^-completeness 
of quantum 3-SAT was recently proven by Nagaj [GN13]. 

This paper is concerned with the quantum 2-SAT problem, which we will also denote simply by 
Q2SAT. One major result concerning this problem is due to Bravyi [Brail], who has proven that it 
belongs to the complexity class P. More precisely, he has proven that Q2SAT can be decided by a 
deterministic algorithm in time O(n^), together with a ground state that has a polynomial classical 
description. In the case of Q2SAT, the Hamiltonian is given as a sum of 2-qubits projectors; each 
projector is defined on a 4-dimensional Hilbert space and can therefore be of rank 1, 2 or 3. In this 
paper, we give an algorithm for Q2SAT of linear complexity. 

Theorem 1 There is a deterministic algorithm for Q2SAT whose running time is 0{n-\-m) where 
n is the number of variables and m is the number of local terms in the Hamiltonian. 

Our algorithm shares the same trial and error approach of the Davis-Putnam procedure for 
classical 2SAT, but handles many difficulties arising in the quantum setting. Eirst, a ground state 
of Q2SAT input may be entangled, some distinctive feature that classical 2SAT does not have. 
Thus the idea of setting some qubit to certain state and propagating from there does not have 
foundation at the first place. Indeed, if a rank-3 projection leaves the only allowed state entangled, 
then any ground state is entangled in those two qubits. We overcome this by showing a product- 
state theorem, which asserts that for any frustration-free Q2SAT instance H that contains only 
rank-1 and rank-2 projectors, there always exists a ground state in the form of a tensor product of 
single-qubit states. 

This structural theorem grants us the following approach: We try some candidate solution 
on a qubit i, and propagate this along the graph. If no contradiction is found, it turns out that 
we can detach the explored part and recurse on the rest of the graph. If a contradiction is found, 
then we can identify two candidates (i, \'tp)i) and (j, \ (f)j) such that either assigning to qubit i 
or assigning \4>)j to qubit j is correct, if there exists a solution at all. More details follow next. 
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To illustrate the main idea of our algorithm, let us suppose that the input contains only pro¬ 
jectors of rank at most two. Such a system can be further simplified to a system consisting only 
of rank-1 projectors, by writing every rank-2 projector as a sum of two rank-1 projectors. Con¬ 
sider, for example, qubits 1 and 2 and a rank-1 projector ni 2 = |V’)(V'I over these qubits. The 
product-state theorem implies that it suffices to search for a product ground state. Thus on the 
first two qubits, we are looking for states |q;), |/3) such that ((a| (8) (/3|)ni2(|a) <8) |/3)) = 0, which 
is equivalent to (a| (8) (/3| • I'*/’) = 0. In other words, we look for a product state |a) (8) |/3) that is 
perpendicular to |^). Assume that we have assigned qubit 1 with the state |a) and we are looking 
for a state |/3) for qubit 2. The crucial point, which enables us to solve Q2SAT efficiently, is that 
just like in the classical case, there are only two possibilities: (i) for any |/3), the state |a) (8) |/3) is 
perpendicular to |'0), or (ii) there is only one state |/3) (up to an overall complex phase), for which 
((a| (8) (/?|) • I'i/’) = 0. The first case happens if and only if |V’) is by itself a product state of the 
form 1^/)) = la"*") (8> |^), where la"*") is perpendicular to |a) and |^) is arbitrary. If the second case 
happens, we say that state |a) is propagated to state |/3) by the constraint state 

The above dichotomy enables us to propagate a product state |s) on part of the system until 
we either reach a contradiction, or find that no further propagation is possible and we are left 
with a smaller Hamiltonian Hg. This smaller Hamiltonian consists of a subset of the original 
projectors, without introducing new projectors. It turns out that once an edge is checked for 
potential propagation, then no matter whether a propagation happens along the edge, the edge 
can be safely removed without changing the satisfiability. Thus the satisfiability of the original 
Hamiltonian H is the same as that of the smaller Hamiltonian Hg. 

We still need to specify how the state |a) is chosen to initialize the propagation. An idea is to 
begin with projectors \'4’){'tp\ for which I'i/;) is a product state |a) (8) |/3). In such cases a product 
state solution must either have !«■*■) at the first qubit or |/?■’■) at the second. To maintain a linear 
running time, we propagate these two choices simultaneously until one of the propagations stops 
without contradiction, in which case the corresponding qubit assignment is made final. If both 
propagations end with contradiction, the input is rejected. 

The more interesting case of the algorithm happens when we have only entangled rank-1 pro¬ 
jectors. What should our initial state be then? We make an arbitrary assignment (say, |0)) to any 
of the still unassigned qubits and propagate this choice. If the propagation ends without contra¬ 
diction, we recurse. If a contradiction is found then we confront a challenging problem. In the 
classical case we could reverse our choice, say xq = 0, and try the other possibility, x* = 1. But in 
the quantum case we have an infinite number of potential assignment choices. The solution is found 
by the following observation: Whenever a contradiction is reached, it can be attributed to a cycle 
of entangled projectors in which the assignment has propagated from qubit i along the cycle and 
returned to it with another value. Then using the techniques of ‘sliding’, which was introduced in 
Ref. [JWZll], one can show that this cycle is equivalent to a system of one double edge and a ‘tail’ 
(see Fig. 1). Using a simple structure lemma, we are guaranteed that at least one of the projectors 
of the double edge can be turned into a product state projector, which, as in the previous stage, 
gives us two possible free choices. 

Let us state here that our algorithm works in the algebraic model of computation: we suppose 
that every arithmetic operation on complex numbers can be done in unit time. There are several 
ways to work in a more realistic model. Bravyi [Brail] suggests considering bounded degree 
algebraic numbers, in which case the length of the representations and the cost of the operations 
should be analyzed. Another possibility would be to consider complex numbers with bounded 
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Figure 1: Handling a contradicting cycle: we slide edges that touch i along the two paths to j 
until we get a double edge with a ‘tail’. We then use a structure lemma to deduce that at least one 
of these edges can be written as a product projector. 


precision in which case exact computation is no more possible and therefore an error analysis 
should be made. By choosing a more convenient computational model, we disregard these issues. 

Classically, Davis-Putnam [DP60] and DPLL algorithms [DLL62] are widely-used heuristics, 
forming the basis of today’s most efficient solvers for general SAT. For quantum k-SAT, it could 
also be a good heuristic if we try to find product-state solutions, and in that respect our algorithm 
makes the first-step exploration. 

Simultaneously and independently from our work and approximately at the same time, de 
Beaudrap and Gharibian [dBG15] have also presented a linear time algorithm for quantum 2SAT. 
The main difference between the two algorithms is how they deal with instances with only entangled 
rank-1 projectors. Contrarily to us, [dBG15] handles these instances by using transfer matrix 
techniques to find discretizing cycles [LMSSIO]. 

2 Preliminaries 

2.1 Notation 

We will use the notation [n] = {1,..., n}. For a graph G = {V, E), and for a subset U C V oi the 
vertices, we denote by G{U) the subgraph induced by U. Our Hilbert space is defined over n qubits, 
and is written as V. = "Hi 01^2 0 • • • 0 'Hn, where T-Li is the two-dimensional Hilbert space of the 
qubit. We shall often write \a)i to emphasize that the 1-qubit state \a) lives in %i. Similarly, \'ip)ij 
denotes a 2-qubit state that lives in Hi 0 Hj. For a 1-qubit state |a) = ao|0) -|- ai|l), we define its 
perpendicular state as la-*-) = ai|0) — ao|l)- 

We shall denote local projectors either by H^, or by Hg, where e = (i, j). When i < j, liij is a 
2-local projector on the qubits i,j; it can be written as Hjj = liij 0 \est^ where Hjj is a projector 
working on Hi 0 Hj and Irest is the identity operator on the rest of the system. Similarly, when 
i = j, liii = tin 0 Irest) where tin is a projector dehned in Hi. Often, in order not to overload 
the notation, we shall use Hjj instead of Hjj, even when acting on states in Hi 0 Hj. Similarly, 
with a slight abuse the notation, we define the rank of a projector Hg to be the dimension of the 
subspace that its local projector Hg projects to, and it will be denoted by rank(ng). We call a 
rank-1 projector Hg = entangled if |'0) is an entangled state, and product if |^/:) is a product 

state. 
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2.2 The Q2SAT problem 

A quantum 2-SAT Hamiltonian on an n-qubit system is a Hermitian operator H = for 

some I C {{i,j) G [n] x [n] : 1 < i < j < n}. We suppose that rank(njj) = 1, for all {i,i) G I, and 
0 < rank(njj) < 4, for all {i,j) G I when i < j. The single-qubit projectors of H as well as its 
2-qubit projectors of rank-3 are called maximal rank. 

The ground energy of a Hamiltonian H = is its smallest eigenvalue, and a ground 

state of H is an eigenvector corresponding to the smallest eigenvalue. The subspace of the ground 
states is called the ground space. A Hamiltonian is frustration-free if it has a ground state that 
is also simultaneously the ground state of all local terms. As explained in the introduction, if the 
Hamiltonian is made of local projectors, it is frustration-free if and only if there is a state that 
is a mutual zero eigenstate of all projectors, which happens if and only if the ground energy is 0. 
Therefore, if |r) is a ground state of a frustration-free quantum 2-SAT Hamiltonian, HelT) = 0 
for all e G I. We can also view each local projector as a constraint on at most two qubits, then a 
ground state satisfies every constraint. 

It turns out that for the representation of the 2QSAT Hamiltonian, it will be helpful to eliminate 
the rank-2 projectors by decomposing each one of them into a sum of two rank-1 projectors. For 
every (i, j) G / such that rank(njj) = 2, let Hjj = Hj^q -|- Hj^-q, where Hjjq and Hj^-q are rank-1 
projectors. Such projectors can be found in constant time. We therefore suppose without loss of 
generality that H is specified by 

H= (hlijq-|-Hjjq), 

rank(nij)7^2 rank{ny)=2 

which we call the rank-1 decomposition of H. 

To the rank-1 decomposition we associate a weighted, directed multigraph with self-loops 
G{H) = (V, E, w), the constraint graph of H. By definition V = {i G [n] : 3j G [n] such that {i,j) ^ 
I or {j,i) G /}, For every rank-3 and rank-1 projector acting on two qubits, there is an edge in 
each direction between the two nodes representing them. For every projector acting on a single 
qubit, there is a self-loop. Finally, for every rank-2 projector, there are two parallel edges in each 
direction between nodes representing its qubits. Because of the parallel edges, E is not a subset of 
V X V. Formally, E = EiU E 2 where 

El = {(i,j) G [n] X [n] : {i,j) ^ ^ and rank(ny) G {1,3}, or {j,i) ^ I and rank(njj) G {1,3}}, 
and 

E 2 = {{i,j,b) G [n] X [n] x [2] : {i,j) G I and rank(njj) = 2, or (j, i) G I and rank(njj) = 2}. 

We say that an edge e G E goes from i to j \i e G {{i, j), {i, j,l), {i, j,2)}. For a projector H 
acting on two qubits, we define its reverse projector by W^'^\a)\ld) = n|/3)|a), and for i < j 
and b G [2], we set Hj* = and Hjib = hlT^. Then for an edge its weight is dehned as 

u’ihj) = bij, and analogously for an edge {i,j,b), we set w{i,j,b) = Hijb- 

Figure 2: (a)The constraint graph for Hamiltonian H = ni 2 -|-ni 4 -|-n 23 -|-n 44 where rank(n 44 ) = 
rank(n 23 ) = 1, rank(ni 2 ) = 2 and rank(ni 4 ) = 3 using its rank-1 decomposition (b) The adjacency 
list representation for the constraint graph G{H) 
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(a) 



(b) 


We will suppose that the input to our problem is the constraint graph G{H) of the Hamiltonian, 
given in the standard adjacency list representation of weighted graphs, naturally modified for 
dealing with the parallel edges as shown in Figure 2. In this representation there is a linked list of 
size at most re containing one element for each vertex, and the element i in this list is also pointing 
towards a linked list containing an element for every edge {i,j) or For an edge this 

element contains j, the projector Yiij and a pointer towards the next element in the list, for an edge 
it also contains the value b. The problem Q2SAT is defined formally as follows. 

Q2SAT 

Input: The constraint graph G{H) of a 2-local Hamiltonian H, given in the adjacency 

list representation. 

Output: A solution if H is frustration free, “Ff is unsatisfiable” if it is not. 


2.3 Simple ground states 

Our algorithm is based crucially on the following product state theorem, which says that any 
frustration-free Q2SAT Hamiltonian has a ground state which is a product state of single qubit 
and two-qubit states, where the latter only appear in the support of rank-3 projectors. A slightly 
weaker claim of that form has already appeared in Theorem 2 of Ref. [000”^ 11]. The difference 
here is that we specifically attribute the 2-qubits states in the product state to rank-3 projectors. 
Just as in Ref. [CCD'^11], our derivation begins with Theorem 1 of Ref. [CCD'^11], which we give 
below. It relies on the notion of a genuinely entangled state in an re-qubit system, which is a pure 
state that is not a product state with respect to any bi-partition of the system. Then Theorem 1 
in [CCD+11] states 

Propositiou 2 A 2-loeal frustration-free Hamiltonian on re qubits which has a genuinely entangled 
ground state always has a produet ground state, whenever n > 3. 

We will also need the following simple fact about 2-dimensional subspaces in (8> 
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Fact 3 Any 2-dimensional subspace V of the 2-qubit space (8) contains at least one product 
state, which can he found in constant time. 

Proof Take a basis {lip), |i?i>)} of the two-dimensional subspace the orthogonal complement of 
V. Our goal is to find a product state |a) 0 |/3) such that (g) |/3)) = (</>|(|a) <8) |/3)) = 0. To 

that aim, expand in the standard basis: \^p) = Ylij'P’ijHd)^ \P^) — Ylij ^ I®) = 

1/3) = Then our task is to find coefficients ai and /3j such that Ylij P^*ij ' — 0 

y',-.-'0*,- • aiBj = 0. We can pass to a matrix notation, in which </>*,• are the entries of 2 x 2 
matrices 'k, <k, and a*, /3j are the coordinates of the 2-vectors a, /3. In that notation, we are looking 
for vectors a, /3 such that 

= 0 . ( 1 ) 

If the matrix <I> is singular, we pick /3 inside its the null space, and choose a such that a^T/3 = 0. 
Otherwise, when $ is non-singular, we let /3 be an eigenvector of the matrix i.e., <k“^'I'/3 = c/3, 

where c is some eigenvalue. Then \k/3 = c<I>/3, and therefore to satisfy Eq. (I), we can choose a such 
that = 0. □ 

Our product state theorem is stated as follows. 

Theorem 4 Any frustration-free Q2SAT Hamiltonian has a ground state which is a tensor product 
of one qubit and two-qubit states, where two-qubit states only appear in the support of rank-3 
projectors. 

Proof Consider a frustration-free Q2SAT Hamiltonian H and let jT) be its ground state. Generally, 
|r) can be written as a product state 

|r) = |q;T)) (g) ® ® , 

where each is a genuinely entangled state defined on a subset 5*-*^ of qubits. Notice if 

Hjj = I — \'p){'ip\jk is a rank-3 projector then necessarily every ground state of H will contain 
\'ip)jk at a tensor product with the rest of the system. Therefore, if happens to be entangled, 
there would necessarily be a subset 5^*^ = {j, k} in the above decomposition with = \'p)jk- On 
the other hand, if \'p)jk is a product state, there would be two subsets and = {k}. 

Let be the Hamiltonian that is the sum of all the projectors whose support is in 

Clearly, is a ground state of H^'^\ By the reasoning in the paragraph above, it is clear that 

for subsets with two or more qubits, that do not correspond to the support of rank-3 projectors, 
the corresponding consists only of rank-1 and rank-2 projectors. It is easy to see that Fact 3 
implies that for such Hamiltonians, which do not contain rank-3 projectors. Proposition 2 also for 
n = 2 case. Therefore, any such also has a ground state 1/3^*)) which is a product state of one 
qubit states: 

|/3«) = |/3«)®|/3«)®... (2) 

The remaining subsets correspond either to one qubit subsets, or to 2-qubits subsets of entan¬ 
gled rank-3 projectors. In all these cases, we define |/3^*^) = 

We now claim that the state |/3) = 1/3^)) (g ... (g) 1/3^)), which is a product of one-qubit and 
two-qubit states, is a ground state of H. To prove this we need to show that this state is in the 
ground space of every projector He in H. If the support of He is inside one of the Si subsets, then 
by definition ne|/3(*^) = 0 and therefore also ne|/3) = 0. Assume then that He is supported on a 
qubit from and a qubit from with i ^ j. We now consider 3 cases: 
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1. If both and contain only one qubit then Ilel/J^*^) 0 = IlelaW) (g) = 0. 

2. If 5(*) is made of one qubit but has two or more qubits, then consider the Schmidt 

decomposition = Ai|xi) 0 \yi) + A 2 IX 2 ) 0 \y 2 )- Here, |xi), \x 2 ) are defined on the qubit 

of Sj that is in the support of He, while |yi), \y 2 ) are defined on the rest of the qubits in Sj. 
The Schmidt coefficients Ai,A 2 are by assumption non-zero, as \aj) is entangled. Then the 
condition Hgla^*^) 0 = 0 is equivalent to Ai|yi) 0 (Hela^*^) 0 kl)) + A 2 |y 2 ) 0 (HelaW) (g) 

\x 2 )) = 0 , and by the linear independence of |yi),|y 2 )) we conclude that Hela^* )) 0 |xi) = 
Hela^*^) 0 \x 2 ) = 0. Therefore, He annihilates the subspace 0 of the two qubits that 
it acts on, and in particular it annihilates 1/3^*^) 0 since |/I«) = |««). 

3. The third case in which both and contain two or more qubits cannot happen. Indeed, 

in such case we write both in their Schmidt decomposition, and from a similar 

argument that was used above, we conclude that He must annihilate 4 independent vectors. 
It therefore cannot be a rank-1 or a rank-2 projector. 

This completes the proof of the theorem. □ 

2.4 Assignments 

Let H = Yle&i ^ 2-local Hamiltonian. By Theorem 4, if H is frustration free then it has a 
ground state which is the tensor product of 1 -qubit and 2 -qubit entangled states, where the latter 
only appear in pairs of qubits corresponding to rank-3 projectors. To build up a ground state of 
such form, our algorithm will use partial assignments, or shortly assignments. An assignment s 
is a mapping from [n]. For every i G [n], the value s(i) is either a 1-qubit state |a), or a 2-qubit 
entangled state \^)ij for some j / i, or a symbol from the set {0, A}. If s{i) = |q;) or s(i) = 
then this value is assigned to qubit variable i, and in the latter case the entangled state is shared 
with variable j. The symbol 0 is used for unassigned variables, and the symbol X is used when 
several values are assigned to some variable. 

We define the support of s by supp(s) = {i G [n] : s{i) 7 ^ 0 }. The assignment s is empty if 
supp(s) = 0. When there is no danger of confusion, we will denote the empty assignment also by 0 . 
We say that an assignment is coherent if for every i, we have s{i) 7 ^ X, and whenever s{i) = 
we also have s{j) = \^)ji. For coherent assignments s and s', we say that s' is an extension of s, 
if for every z, such that s{i) 7 ^ 0 , we have s'{i) = s{i). A coherent assignment is total if s{i) 7 ^ 0 , 
for all i. Clearly, a coherent assignments defines a product state of 1-qubit and 2-qubits states on 
qubits in its support. We denote this state by |s). We say that a coherent assignment s satisfies 
a projector He, or simply that it satisfies the edge e, if for any total extension s' of s we have 
Hels') = 0 . 

For H = X^es/given in rank-1 decomposition, and a coherent s, we define the reduced 
Hamiltonian Hg of s as 

Hg = H- ^e. 

s satisfies e 

We will denote the constraint graph G{Hs) of the reduced Hamiltonian Hg by Gg = (VgjEg). We 
call a coherent assignment s a pre-solution if it has a total extension s' satisfying every constraint in 
H, and we call s a solution if s itself satisfies every constraint in H. Obviously, an assignment is a 
solution if and only if Gg is the empty graph. A coherent assignment s is closed if supp(s) HI 4 = 0. 


3 Propagation 


The crucial building block of our algorithm is the propagation of values by rank-1 projectors. This 
is the quantum analog of the classical propagation process when for example the clause Xi V Xj 
propagates the value Xj = 0 to the value Xj = 1 in the sense that given Xj = 0 , the choice Xj = 1 is 
the only possibility to make the clause true. In the quantum case this notion has already appeared 
in Ref. [LMSSIO], and can in fact be traced back also to Bravyi’s original work. Here, we shall 
adopt the following definition 

Definition 1 (Propagation) Let Hg = |V')(V'| a rank-1 projector acting on variables i,j, and 
let la) be either a l-qubit state assigned to variable i, or a 2-qubit entangled state assigned to 
variables k, i for some k ^ j. We say that Hg propagates |a) if, up to a phase, there exists a unique 
l-qubit state \f3) such that Hgja) 0 \l3)j = 0. In such case we say that |a) is propagated to |/3) along 
Hg, or that Hg propagated |a) to \j3). 

The following lemma shows how the propagation properties of Hg = |?/>)('!/’| are determined by 
entanglement in \'ip). 

Lemma 5 Consider the rank-1 projector He = defined on qubits i,j. If {ip) is entangled, it 

propagates every l-qubit state \a)i to a state \(3{a))j such that if\a)i / |a')i then |/3(a))j 7 ^ |/3(a'))j. 
This propagation can be calculated in constant time. When \fj) is a product state \ip) = |x)j ( 8 ) \y)j, 
the projector Hg does not propagate states that are proportional to Ix"*-)*, while all other states are 
propagated to \y'^)j. 

Proof Assume that \'tp) is entangled and consider the state |a). Our task is to show that there 
always exists a unique |/3) (up to an overall constant) such that n(|a) ( 8 ) |/3)) = 0 , and that different 
|a)’s yield different \(3)’s. 

Expanding {fj), |a), and |/3) in the standard basis {ip) = ® b')! 1*^) = 

1/3) = Pjlj), the condition Hgda) (8) |/3)) = 0 translates to Yij — 0- Assuming that \'ip) 

is entangled, one can easily verify that the 2x2 matrix (V’))) non-singular. Then using the simple 
fact that in a two-dimensional space every non-zero vector has exactly one non-zero vector (up to 
an overall scaling) which it is orthogonal to, it is straightforward to deduce that for every non-zero 
vector (ao, ai) there is a unique (up to scaling) non-zero vector (/3o, /3i) such that Yi j ~ 0- 

Moreover, (/3o,/3i) can be calculated in constant time, and that different (ao,ai) necessarily yield 
different (/3o,/3i). 

The case when \ip) is a product state is straightforward. □ 

We now present two lemmas that describe the structure of the global ground state of the 
system, if we know that part of it is in a tensor product of l-qubit or 2 -qubits states, which are 
then propagated by some Hg. 

Lemma 6 (Single qubit propagation) Consider a frustration-free Q,2SAT system H = Z^gg/Hg 
with a rank-1 projector Hg = \'ip){'tp\ between qubits i,j, and assume that H has a ground state of 
the form jT) = \a)i 0 \rest). Then: 

1- //Hg propagates \a)i to |/3)j then necessarily \rest) = \ f3)j ® \rest'). 

2. |r) is also a ground state of the Q2SAT Hamiltonian H — He. 
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Proof For the first claim assume that He propagates \a)i to |/?)j. Without loss of generality, we 
may expand 

\rest) = |/3)j (g) \resti) + \j3^)j ® \rest 2 ) , 

where the states |resti), |rest 2 ) are defined on all the qubits of the system except for (i, j), and are 
not necessarily normalized. Plugging this expansion into the condition Heir) = 0, we obtain the 
equation 

{I\.e\a)i\l3)j) ® \resti) + {I\.e\a)i\l3^)j) ® \rest 2 ) = 0. 

Since He propagates \a)i to \(3)j, we have He\a)i\l5)j = 0 and He \a)i\(3-^)j / 0. Therefore, the 
above equation implies that |rest 2 ) = 0 , and we may set \rest') = \resti). 

The second claim follows trivially from the frustration-freeness of the system. □ 

Lemma 7 (Entangled 2-qubits propagation) Consider a frustration-free Q2SAT system H 
with a rank-1 projector He = |'0)(V’l between qubits i,j. Assume that H has a ground state of 
the form |r) = \4>)ik <Si \rest), where |(/>) is an entangled state on qubits i, k with k ^ j. Then: 

1. \'4:) is a product state IV') = |a^)|y)- 

2. He propagates \4>) to \y-^) and necessarily \rest) = \y'^)j ® \rest'). 

3. |r) is also a ground state of the Q2SAT Hamiltonian H — He- 

Proof Write |0) in its Schmidt decomposition |^) = Ai|a) ® |/3) + A 2 |a"'') < 8 ) and note that 
both Ai, A 2 / 0, since |</>) is entangled. Plugging this into the condition flelP) = 0, we get 

rieir) = Xi\/3)k < 8 ) ne(|a)j (g |rest)) + X2\/3-^)k ® ne(|Q:-^)i ( 8 ) \rest)) = 0 . 

Since \(5) is is linearly independent of l/S-*-), we conclude that ne(|a)i( 8 )|rest)) = ne(|a“'“)i( 8 )|rest)) = 

0 . 

To prove the hrst claim, assume by contradiction, that If:) is entangled. Then by Lemma 5, lie 
propagates |a) and la-*-) to two different states, say, I 71 ) / I 72 ). But then by Lemma 6 , it follows 
that I rest) must be both in the form | 7 i)j ( 8 ) |rest') and | 72 )j ® |rest') - which is a contradiction! 

For the second claim, assume that lip) = jx) ( 8 ) \y) is a product state. Since ne(|a)i ( 8 ) |rest)) = 
He (|a-*-)i ( 8 ) |rest)) = 0 , both states \a)i ® |rest), ( 8 ) |rest) are ground states of the single 
projector Hamiltonian H = Hg. Using Lemma 5 and Lemma 6 , together with the fact that and at 
least one of the states la), la-*-) is different from |a:"''), we conclude that |rest) = |y'')j ( 8 ) |rest'). 
The third claim, as before, follows simply from the frustration-freeness of the system. □ 

Let FI be a Q2SAT Hamiltonian in rank-1 decomposition, let s be a coherent assignment, and let 
Gg = {Vs, Eg) be the constraint graph of the reduced Hamiltonian Hg. We would like to describe in 
Gg the result of the iterated propagation process when a value given to variable i is propagated along 
all possible projectors, then the propagated values are propagated on their turn, and so on until 
no more value assigned during this process can be further propagated. The propagation can get 
started when the initial value is already assigned by s, that is when s(i) = |(5) for |(5) G {|a), | 7 )ij}, 
where |a) is some 1 -qubit state and \^)ij some a 2 -qubit state, or it can get started when s{i) = 0, 
in which case we shall explicitly choose a 1 -qubit state |a) and assign it to f. 
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Let now s,i and |(5) be such that s{i) G {0, |(5)}. We say that in the constraint graph Gg an 
edge e G Eg from i to j propagates |(5) if lie propagates it, and we denote by prop(s, e, |(J)) the 
state Id) is propagated to. We generalize the notion of propagation in Gg from edges to paths. Let 
i = io 0 i, • • • be vertices in 14, and let ej be an edge from ij to ij+i, for j = 0,..., /c — 1. Let 
s{i) G {0, |(5)}, and set |ao) = li^)- Let |ai),..., \ak) be states such that the propagation of \aj) 
along Ile^ is |aj_|_i), for j = 0,..., /c — 1. Then we say that the path p = (eo,..., ek-i) from zq to 
ik propagates |(5), and we set prop(s,p, |(5)) = \ak)- We say that a vertex j G 14 is aecessible by 
propagating 1 5) from i if either j = i or there is a path from i to j that propagates 1 5). We denote 
by 1 5)) the set of such vertices, and by exts'^°^(z, |(5)) the extension of s by the values given 

to the vertices in |5)) by iterated propagation. 

Let us suppose that s' = exts'’°^(i, |(5)) is also coherent. The set |(5)) divides the edges 

Eg into three disjoint subsets: the edges Ei of the induced subgraph |(5))), the edges 

E 2 between the induced subgraphs G{yr°'‘^{'i, |(5))) and G{Vg \ |(5))), and the edges Lis of 

the induced subgraph G{Vg \ Id))). While the edges in Ei U E 2 are satisfied by s', none of 

the edges in E^, is satisfied. Therefore Gg/ is nothing but GiVg \ |5))) without the isolated 

vertices, and it can be constructed by the following process. Given s and i, the edges in E 1 UE 2 can 
be traversed via a breadth first search rooted at i. The levels of the tree are decided dynamically: 
at any level the next level is composed of those vertices whose value is propagated from the current 
level. The leaves of the tree are vertices in 14 \ |(5)). The algorithm Propagation uses a 

temporary queue Q to implement this process. 


Algorithm 1 Propagation(s, G.,41-5)) 
s(z) := |5) 

create a list L and a queue Q, and put i into Q 
while s is coherent and Q is not empty do 
remove the head j of Q 
for all edge e from j to k do 
remove e from Eg 
if e propagates s{j) then 

prop(s,e,s(j)) if 
X if 

enqueue k 

if e is not propagating and k is 
remove j from I 4 

if s is not coherent return “unsuccessful” 
for all k in L do 

for all edges e from k to i do 

if £ was removed from I 4 then remove e 
if all edges outgoing from k were removed then remove k from I 4 



s{k) = 0 

s{k) 0 {0,prop(s,e,s(j))} 
not in L then put k into L 


Lemma 8 (Propagation Lemma) 

Let Propagation(s, Gs, i, |5)) be ealled when Hg doesn’t have rank-3 constraints, and s{i) G 
{ 0 , |(5)}. Let s' and G' = {V', E') be the outcome of the procedure. Then: 
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1. //Propagation(s, Gs, i, \5)) doesn’t return “unsuccessful” then s' = exts'^°’’(z, |5)) and G' = Gg/. 
Moreover, if s is a pre-solution then s' is a pre-solution, and if s is closed then s' is also closed. 

2. //Propagation(s, Gs, i) returns “unsuccessful” then there is no solution z of which is an ex¬ 
tension of s and for which z{i) = |(5). 

3. The complexity of the procedure is 0{\Es\ — \Es'\). 


Proof The assignments made during the breadth first search correspond exactly to the the 
paths propagating |(i) from i, therefore the extension of s created by the process is indeed s' = 
exts™^(i, |(5)). The while loop removes the edges between vertices in |(5)) and the edges 

which go from |(5)) to 14 \ Vr°^{i, |(5)), as well as the vertices in Vr°^{i-, |5))- Then the 

edges from 14 \ 14'’'^°'’(4 l<5)) to |5)) are removed, as well as the remaining vertices without 

outgoing (and incoming) edges. Therefore we have G' = Gg'- 

Let us suppose that s is a pre-solution, and let z be an extension of s which is a solution and 
which is a product state on the vertices in 14- By Theorem 4 there exists such a solution since Hg 
doesn’t have rank-3 constraints. We define the assignment z' by 


Aj) 


s'{j) if j G supp(s') 
z{j) otherwise. 


Then z' is a solution which is an extension of s', and therefore s' is a pre-solution. If s is closed 
then so is s' since only the vertices in Vr°'^{i, |(5)) get assigned during the process, and they are not 
included into 14 '. 

Let us now suppose that the procedure returns “unsuccessful”. Then there is a vertex k G 
Vr°'^{i, |5)), and two paths p and p' in Gg from itok such that prop(s,p, |(5)) = |/3), prop(s,p', |<5)) = 
1/3') and 1,0) / |/3'). Let us also suppose that there exists a solution z which is an extension of s 
and for which z{i) = |(3). Then by the repeated use of Lemma 6, and also by using once Lemma 7 
when |(5) is a 2-qubit entangled state, we conclude that z{k) is simultaneously equal to |/3) and to 
1/3'), which is a contradiction. 

Finally Statement 3 follows since every step of the procedure can be naturally charged to an 
edge in Eg \ Egi, and every edge is charge only a constant times. 

□ 


4 The main algorithm 

4.1 Description of the algorithm 

We now give in broad lines a description of our algorithm we call Q2SATSolver. It takes as input the 
the adjacency list representation of the constraint graph G{H) of a 2-local Hamiltonian H in rank-1 
decomposition. The algorithm uses four global variables: assignments sq si initialized to 0, and 
graphs Go and Gi in the adjacency list representation, initialized to G{H). The algorithm consists 
of four phases, and except the first one, each phase consists of several stages, where essentially 
one stage corresponds to one Propagation process. In the case of an unsatisfiable Hamiltonian 
the algorithm at some point outputs “H is unsatisfiable” and stops. This happens when either 
the maximal rank constraints are already unsatisfiable, or at some later point several values are 
assigned to the same variable during a necessary propagation process. 
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In the case of a frustration-free Hamiltonian, at the beginning and end of each stage, we will 
have So = si, fand Gq = Gi = Gsq- In the first two phases only (so, Go) develops, and is copied to 
(si, Gi) at the end of the phase. In the last two phases, (sq, Go) and (si, Gi) develop independently, 
but only the result of one of the two processes is retained and is copied into the other variable at 
the end of the phase. This parallel development of the two processes is necessary for complexity 
considerations, it ensures that the useless work done is proportional to the useful work. 

In the first phase the procedure MaxRankRemoval satisfies, if this is possible, all constraints of 
maximal rank. In the second phase all these assignments are propagated, which, if successful, results 
in a closed assignment s such that Hg has only rank-1 constraints. In the third phase the procedure 
ParallelPropagation satisfies the product constraints one by one and propagates the assigned values. 
To satisfy a product constraint, the only two possible choices are tried and propagated in parallel. In 
the fourth phase the remaining entangled constraints are taken care of, again, one by one. To satisfy 
a constraint, an arbitrary value is tried and propagated. In case of an unsuccessful propagation 
we are able to efficiently find a product constraint implied by the entangled constraints considered 
during the propagation, and therefore it becomes possible to proceed as in phase three. In case of 
success we are left with a satisfying assignment and the empty constraint graph. Theorem 1 is an 
immediate consequence of the following result. 


Algorithm 2 Q2SATSolver(G(//)) 

So = si := 0, Go = Gi := G{H) > Initialize global variables 

MaxRankRemoval() i> Remove maximal rank constraints 

while there exist i € Vq such that s{i) ^ 0 do t> Propagate all assigned values 

Propagate(so, Go, f, so{i)) for some vertex i in Go such that so{i) ^ 0 
if the propagation returns “unsuccessful” output is unsatisfiable” 

Si •= So, Gi := Go 

while there exists in Gq a product edge with constraint Ioq )*o 0 I®!")*!(“o'lio 0 (“I'lu 

ParallelPropagation(io, |ao),fi, |cii)) Remove product constraints 

while Go is not empty do > Remove entangled constraints 

ProbePropagation(f) for some vertex i 

output |s) for any total extension s of so. 


Theorem 9 Let G{H) = (P, E) he the constraint graph of a 2-local Hamiltonian. Then: 

1. If H is frustration-free, the algorithm Q2SATSolver(G(i/)) outputs a ground state |s). 

2. If H is not frustration-free, the algorithm Q2SATSolver(G(i?)) outputs “H is unsatisfiable”. 

3. The running time of the algorithm is 0(|P| -|- |ili|). 

Theorem 9 will be proven in Section 4.5. 
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4.2 Max rank removal 

The MaxRankRemoval procedure is conceptually very simple. Since every maximal rank constraint 
has a unique solution (up to a global phase), it makes this assignment for each constraint, and then 
checks if this is globally consistent. 


Algorithm 3 MaxRankRemoval() 


for all i G Vo such that rank(njj) = 1 and \4>) is the unique state satisfying Ilj. 
so(i) := |(/)) 

for all i £ Vq, for all edge e G Eq from i to j such that rank (lie) = 3 do 
if \a)i\f3)j is the unique product 2-qubit state satisfying He then 

\\a) ifso(*) = © 

■SnU) := \ 

\X if so(i)^ {©,!«)} 

if 1 7) is the unique entangled 2-qubit state satisfying He then 
\-i)ij if so(*) = © 

X if so(*) i {©, 


so{i) := 


if So is not coherent for some z G 14 or for three distinct variables we have s(i) = \'y)ik 

and He is an entangled projector on i and j then 
ontput “41 is unsatisfiable” 

remove from Eq every edge e such that He is satisfied by sq- 
remove every isolated vertex from Go 
Si •= So, Gi := Go 


Lemma 10 Let so,Gq,si,Gi be the outcome 0 /MaxRankRemoval. Then: 

1. If MaxRankRemoval doesn’t output “H is unsatisfiable” then sq is coherent, it satisfies every 
maximal rank constraint, Go = G{Hsg) and so = si,Go = Gi. Moreover, if H is satisfiable 
then So is a pre-solution. 

2. If MaxRankRemoval outputs “H is unsatisfiable” then indeed H is unsatisfiable. 

3. The complexity of the procedure is 0(|1/| -|- |41|)|. 

Proof If the procedure doesn’t output “44 is unsatisfiable” then indeed so is coherent and it 
satisfies all maximal rank constraints. The removal of the necessary edges and vertices insures that 
Go = G(44 sq), and obviously so = si,Go = Gi. If 44 is satisfiable, then it has a ground state for 
some total assignment s. This s is an extension of so because there is a unique way to satisfy the 
maximal rank constraints. 

Maximal rank projectors are such that there is a unique assignment for their qubit(s) which 
satisfies them. The first part of the procedure creates the assignment which assigns these necessary 
values. If this assignments is not coherent then H is unsatisfiable. Similarly, if so assigns an 
entangled 2-qubit state between variables i and k, and there is an entangled rank-1 constraint 
between i and j, then by Lemma 7 it is impossible to extend so into a satisfying assignment, and 
therefore 44 is unsatisfiable. This proves Statement 2. 
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The procedure can be executed by a constant number of vertex and edge traversals for sq, and 
similarly for si. □ 

4.3 Algorithm ParallelPropagation 

The procedure ParallelPropagation is called when sq is a closed assignment, and in Gsq there is a 
product edge. Since there are only two ways to satisfy a product constraint, these are tried and 
propagated in parallel. If one of these propagations terminates successfully, the other is stopped, 
which ensures that the overall work done is proportional to the progress made. 


Algorithm 4 ParallelPropagation(io, |ao), h, |ai)) 

Run in parallel Propagation(so, Go, |cro)) and Propagation(si, Gi, zi, |cri)) 
until one of them terminates successfully or both terminate unsuccessfully 

if both propagations terminate unsuccessfully then 
ontput “R is unsatisfiable” 

else let Propagation(so, Go, zo, |<ao)) terminate first (the other case is symmetric) 
undo Propagation(si, Gi, zi, |ai)) 

Si := So, Gi := Go 


Lemma 11 Let ParallelPropagation be called when sq is closed, doesn’t have rank-3 constraints, 
Go = Gsq, in Go there exists a product edge from zq to i\ with constraint [oq ) 0 loi"), and si = 
So, Gi = Go. Let Sq, s)^, Gg, G'l be the outcome of the procedure. Then: 

1. //ParallelPropagation doesn’t output “H is unsatisfiable” then Sg is a proper closed extension 
of So, Gg = Gg'^, and s'^ = Sg, G'^ = Gg. Moreover, if s is a pre-solution then Sg is a 
pre-solution. 

2. If ParallelPropagation outputs “H is unsatisfiable” then indeed H is unsatisfiable. 

3. The eomplexity of the procedure is O(|//so| — 

Proof If the procedure doesn’t output “// is unsatisfiable” then at least one of the parallel 
propagations terminates successfully, say Propagation(so, Go, zo, |«o))- Then Sg is a proper extension 
of So since sq is closed and therefore so(zo) = ©• Obviously s'^ = Sg and G'l = Gg, and all other 
claims follow from the Propagation Lemma. 

Since HgQ doesn’t have rank-3 constraints, by Theorem 4 if it is frustration free, it has a product 
ground state. In Hg^ there exists a product edge from zq to zi with constraint |ag )io <8> Icri")*! (cio U ® 
(«;!'■ Itherefore only assignments which have either |ao) assigned to variable zo or |ai) assigned 
to variable ii can be a solution. But if both propagations output “unsuccessful”, then by the 
Propagation Lemma no such assignment can satisfy Hg^. Therefore Hg^ is not frustration free, and 
neither is H. 

For the complexity analysis observe that the unsuccessful or unterminated propagation of the 
parallel processes makes at most as many steps as the successful one. This is the reason for 
performing the two propagations in parallel. Undoing this propagation can be performed in the same 
order of time as the propagation itself, for example, by copying the removed edges into temporary 


15 






Figure 3: The sliding of the edge (io)*i) over the path —)■ ^2 —... —> ifc, until it becomes the 

edge (io,4)- 


lists. The claim on the complexity of the successful propagation follows from the Propagation 
Lemma. □ 

4.4 Algorithm ProbePropagation 

The procedure ProbePropagation is evoked when sq is a closed assignment, and in Gs^ there are 
only entangled constraints. It picks an arbitrary vertex va. i G Vg, assigns |0) (an arbitrary value) 
to it, and propagates this choice. In the lucky case of successful propagation this is repeated. 
Otherwise, we reach a contradiction: there is some j G 14, such that two propagating paths assign 
different values to it. We prove below the Sliding Lemma which already appeared in Ref. [JWZIl]. 
It implies that when zq 4 is a propagating path of entangled rank-1 projectors, the 

ground space of the Hamiltonian + • ■ ■ + i® equal to the ground state of the 

Hamiltonian + • • • + where is a new projector defined on the qubits 

(io,4) that replaces the projector Graphically, this can be viewed as if we are sliding the 

(io, ii) edge on the path ii as shown in Fig. 3. Therefore, if we have two propagating 

paths starting at i and ending at j, they define two projectors on qubits as illustrated in 

Fig. 1. As we shall see, if these two paths are contradicting then necessarily the two projectors are 
different, which by Lemma 3 implies the existence of a product constraint between {io,ik) variables. 
In such case, we can proceed by calling the procedure ProbePropagation. 

Lemma 12 (Sliding Lemma) Consider a system on 3 qubits i,j and k. Suppose that we have 
a two rank-1 constraints Hi = {tpi){ipilij on qubits {i,j) n 2 = \tp 2 ){'ip 2 \jk on qubits {j,k). If 
|?/: 2 ) is entangled, there is another rank-1 constraint Ha = \'if 3 ){'tp 3 \ik on qubits {i,k) such that the 
ground space of Hi -|- n 2 is identical to the ground space ofIl 2 + Ha. In addition, if a single qubit 
state \a)i is propagated by Hi -|- n 2 to \l3)k, then it is also propagated to \/3)k directly via Ha. 

Proof Consider the Schmidt decomposition \if 2 )jk = ^i\xi)j\yi)k + ^ 2 \x 2 )j\y 2 )k-, where Ai, A 2 / 0, 
as \'4^2)jk is entangled. Define a non-singular transformation T on qubit j by AiTlxi) = \y 2 ) and 
A 2 T|® 2 ) = -\yi)- Then T\'ip2)jk = \y2)j\yi)k - \yi)j\y2)k is the anti-symmetric state. Let \'^i)ij 
and \'if 2 )jk be the normalization of r|V’i)*j and T\'il) 2 )jk respectively, and use them to define the 
rank-1 projectors Hi, n 2 . Since n 2 projects into the anti-symmetric subspace, then any state in the 
ground space of Hi -|- n 2 must be invariant under a swapping of qubits j, k. Therefore, definining 
\'4’3)ik = and Ha = I — |'03)ifc('03|*fc, the the ground space of Hi -|- n 2 is identical to the 

ground space of Ha -|- Ll 2 - Applying now the inverse transformation T~^ on qubit j, the projector 
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n2 returns to 112, while IIs remains unchanged. Since both T,T~^ are non-singular, it follows that 
ground space of Hi -|- 112 is identical to the ground space to 112 - 1 - 113 . 

For the second claim, assume by contradiction that Ha does not propagate la)* to \f3)k- Then 
there is a 1-qubit state I 7 ) 7^ |/3), such that n3(|Q!)i|7)fc) = 0. Since 112 is a rank -1 entangled 
projector, by Lemma 5, it propagates \'y)k to some state |(j)j, and therefore the state \c()i\S)j\j)k is 
a ground state of 112 + 1 ^ 3 , as well as of Hi - 1 - 112 . But this contradicts the assumption that latter 
propagates \a)i to \f3)k- D 

Applying Lemma 12 iteratively, we reach the following corollary 

Corollary 13 Let H = X^es/ ® 2-local Hamiltonian in rank-1 decomposition. Let io, ii,... ik 

be vertices in V, and let ej be an edge from ij to ij+i, for j = 0,..., k — 1 such that the rank-1 con¬ 
straints Ylf,. are entangled. Then there exists a 2-qubit entangled state I 7 ) between io and ik such that 
the ground space 0 /X]j=o is identical to the ground space o/X]j=i + l7)(7|io,4- Moreover, 
if is propagated to |/3)ij, along the path, then it is also propagated directly by I 7 ) (71*0,4- 

We will denote the state I 7 ) in the conclusion of the corollary by slide(p). 


Algorithm 5 ProbePropagation(i) 

Propagation(so,Go,i, |0)). 

if the propagation is successful then si := sq, Gi := Gq 
else 

Let j such that |so(j)| > 1 

find two paths pi and p 2 in Gq from i to j such that prop(so,pi, |0)) 7 ^ prop(so,P 2 , |0)) 
find a product state la"*")* 0 \fd'^)j in the two dimensional space span{slide(pi), slide(p 2 )} 
undo Propagation(so, Go, i, |0)) 

ParallelPropagation(i, \a),j, |/3)) 


Lemma 14 Let ProbePropagation be called when so is closed, has only rank-1 entangled con¬ 
straints, Go = Gsq, and si = so,Gi = Go- Let Sg,s'^,Gg,G( be the outcome of the procedure. 
Then: 

1. //ProbePropagation doesn’t output “H is unsatisfiable” then Sg is a proper closed extension of 
So, Gg = Gg'^, and s( = Sg, G( = Gg. Moreover, if s is a pre-solution then Sg is a pre-solution. 

2. If ParallelPropagation outputs “H is unsatisfiable” then indeed H is unsatisfiable. 

3. The complexity of the procedure is O(|//so| — 

Proof If the procedure doesn’t output “// is unsatisfiable” then either Propagation(so, Go, i, |0)) 
or one of the parallel propagations (say Propagation(so, Go, i, |a))) terminates successfully. Then 
Sg is a proper extension of sq since sq is closed and therefore so(io) = ©• Obviously s'^ = Sg and 
G( = Gg, and all other claims follow from the Propagation Lemma. 

Let’s suppose that all three propagations are unsuccessful. By Corollary 13, any solution for 
HgQ also satisfies let"*")* 0 |/3■’■)j(a■’■|^ 0 {/3'^\j. Then Lemma 11 implies that H is unsatisfiable. 

For the complexity analysis the interesting case is when the first propagation, that we call 
Propfaiiure, is unsuccessful but one of the two parallel propagations is successful. Let’s call this 
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successful one Propsuccess- The main observation here is that every propagating edge in Propfaiiure 
will also be propagating in Propsuccess, since by Lemma 5 entangled edges always propagate. The 
paths Pi and p 2 can be found in time proportional to the size of the subgraph visited by Propfaiiure- 
Indeed, observe that the edges of the two paths, except the last edge of one of the two, are edges 
in the BFS tree underlying PropfaUure- The way from a vertex to the root of the tree can be then 
found, for example, by maintaining for each vertex in the tree, a pointer towrads its father. The 
product state |a) 0 |/3) can be found in constant time by Lemma 3. Therefore, by Lemma 11, the 
complexity is indeed 0{\Esq \ — □ 

4.5 Analysis of the algorithm 

Proof of Theorem 9 If Lf is frustration free then by Lemma 10 MaxRankRemoval outputs a 
pre-solution sq that satisfies every maximal rank constraint. By the Propagation Lemma, at the 
end of Phase two, in addition sq is a closed. By Lemma 11 ParallelPropagation outputs sq such that 
in addition in Hg there are only entangled constraints. By Lemma 14 at the end of the algorithm 
in addition Hg is empty, and therefore s is a solution. 

If the algorithm doesn’t output “H is unsatishable” then by Lemma 10, by the Propagation 
Lemma, and by Lemmas 11 and 14 it outputs a coherent assignment s such that Gg is the empty 
graph, and therefore s is a solution. 

The complexity of MaxRankRemoval by Lemma 10 is 0{\E\). After the second phase, the prop¬ 
agation of the assigned values during MaxRankRemoval, the copying of sq and Go into respectively 
Si and Gi can be done by executing the same propagation steps this time with si and Gi. The 
complexity of the rest of the algorithm by the Propagation Lemma, and Lemmas 11 and 14 is a 
telescopic sum which sums up to also 0{\E\). 
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